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1.0 GENERAL INFORMATION 

1 . 1 Summary 

This report provides user information for program PROFILE, 
an aerodynamics design utility for refining, plotting, and 
tabulating airfoil profiles. The theory and implementation 
details for two of the more complex options are also presented. 
These are the REFINE option, for smoothing curvature in selected 
regions while retaining or seeking some specified thickness 
ratio, and the OPTIMIZE option, which seeks a specified 
curvature distribution. REFINE uses linear techniques to 
manipulate ordinates via the central difference approximation to 
second derivatives, while OPTIMIZE works directly with curvature 
using nonlinear least squares methods. 

Use of programs QPLOT and BPLOT is also described, since 
all of the plots provided by PROFILE (airfoil coordinates, 
curvature distributions) are achieved via the general-purpose 
©PLOT utility. BPLOT illustrates (again, via QPLOT ) the shape 
functions used by two of PROFILE'S options. These three 
utilities should be distributed as one package. All are 
implemented in FORTRAN 77. They were developed on a VAX-1 1/780 
under VMS for the Applied Aerodynamics Branch at NASA Ames 
Research Center, Moffett Field, California, by staff from the 
Palo Alto office of Informatics General Corporation, under 
Contract NAS2-11555-307-1 . 

The plots produced by QPLOT require the widely-available 
DISSPLA* graphics package, which may inhibit some potential 
users but PROFILE and BPLOT are readily portable. Translation 
of the plotting part of QPLOT to some other graphics package 
should be straightforward for those without access to DISSPLA. 
Plots are crucial to the effectiveness of PROFILE and BPLOT. 


1.2 Introduction 

Program PROFILE is a utility for computational aerodynamic 
applications. Originally intended for plotting airfoil 
coordinates and tabulating derivatives, it has evolved into a 
more powerful package capable of modifying airfoils in a number 
of ways. It has thus become a useful design tool. Profiles may 
be perturbed through interactive application of shape functions, 
refined by seeking desired second derivative distributions while 

* DISSPLA is a proprietary product of ISSCO, San Diego, and 
is not available from NASA. It may be purchased under 
license . 
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adjusting thickness ratio, or optimized by seeking target 
curvature distributions directly. 


Host of the options apply to a single two-dimensional 
airfoil section per run, although handling of multiple sections 
in the same input data file is provided for in two cases_ 
sections of a wing (new plot page for each section); and the 
elements of a multi-element airfoil such as a wing with a slat 
(all elements appear on the same plot frame) . The 
sophisticated options offer "before" and "after plots of both 
the airfoil and its curvature distributions (upper and lower 
surfaces on separate frames). 


The simpler options include redistributing the data 
in a variety of ways, and renormalizing the coordinates, 
standardized geometry data formats are provided 
compatibility with most computational flow codes. 


points 

Three 

for 


PROFILE provides a cheap and easily-used means of 
illustrating and refining airfoil geometry data. It should 
serve as a flexible vehicle for additional capabilities, vhioh 
airfoil designers are encouraged to incorporate as their needs 

arise . 


1 . 3 Environment 

PROFILE is an interactive program developed on a Digital 
Equipment Corporation VAX-1 1/780 under VMS V3.6. It is written 
in FORTRAN 77, with VAX extensions used in certain places where 
adhering to the strict standard was either out of the question 
(such as the READONLY keyword on explicit OPEN statements for 
reading files not owned by the user), or considered too 
restrictive (as in the case of the long variable names in 
certain routines). All computations are performed in single 
precision (32-bit words on the VAX). 

Plotted results are essential to effective use of PROFILE. 
However plotting is invariably installation-specific. This is 
one reason that PROFILE itself does not perform any graphics - 
it is more portable that way. The other reason is that some of 
the plot table data output from one run of PROFILE (specifically, 
the curvature data) may be reused as input for another run + 
Instead advantage has been taken of a general-purpose plot 
package known as QPLOT , developed under the same contract as 
PROFILE. QPLOT uses DISSPLA (referenced above), and should 
accompany any distribution of PROFILE. 


1 1 ! 
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The "device-independent " intermediate files produced by 
©PLOT (for previewing on a video display terminal, or post- 
processing on a hard-copy device such as a Versatec plotter) are 
in a format peculiar to the NASA Ames environment, although a 
simple change in the device initialization can produce other 
types of intermediate files. See Appendix B for more on QPLOT. 


A second supplementary program that should accompany 
distributed copies of PROFILE is known as BPLOT , described in 
Appendix A. BPLOT produces plots - again, via QPLOT - of the 
"bump" functions employed by PROFILE in its MODIFY and OPTIMIZE 

modes . 


1.4 Future Directions 

An enhanced PROFILE, contemplated at the time of writing, 
would provide for truly interactive graphics as opposed to 
end-of— run graphics. For instance, in the MODIFY mode, 
perturbations could be tried, viewed, and perhaps undone and 
retried, in a single session instead of requiring several runs 
of PROFILE as at present. The graphics could be fairly crude 
(given the resolution of today's video display terminals), would 
probably provide windowing (i.e., enlarging, for the same 
reason), and could be achieved with a much less bulky graphics 
package than DISSPLA. 

Generation of Cp distributions is another possibility 
appropriate to PROFILE given that this can be done cheaply with 
good approximations, at least for subsonic flows. Boundary 
layer displacement thickness distributions could be generated 
similarly . 

Treating curvature as a single distribution rather than 
separate distributions for upper and lower surface may be more 
effective for some purposes. The authors and sponsors of 
PROFILE encourage users to adapt this tool to suit their needs, 
in the belief that a sound framework has been provided in the 
present code for capabilities unforeseen so far. Feedback is 
welcomed, and conscientious adherence to the standards of modern 
software engineering is urged upon those designers inclined to 
perpetuate the programming style (or lack thereof) for which 
engineers have long been notorious. 
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2.0 BASIC INPUT DESCRIPTION 

PROFILE has numerous modes of operation, some of which 
require considerable activity at the keyboard during program 
execution. Details of such interactive entries are described in 
the appropriate subsections of Chapter 4. The basic inputs 
common to all modes are described here. 

An airfoil dataset and control inputs are required for 
every run of PROFILE, in two separate files. PROFILE prompts 
for these file names. 


2.1 Airfoil Data Format 

The airfoil dataset may be input in one of three formats, 
and also may be output in any of the three formats described 
below. The first format shown is the most oommon and is known 
as "STANDARD" or "PROFILE" format. 


LEGEND 

NU UPPER SURFACE 
X Y 

X Y 


NL LOWER SURFACE 
X Y 

X Y 


< CHARACTER* 80 > 

<14, IX, UPPER SURFACE > 

< 2F10 . 0 > 

< 2F10 . 0 > 

(May be X/C, Y/C; 
Xs are increasing . ) 

<14, IX, LOWER SURFACE > 

< 2F10 . 0 > 

< 2F10 . 0 > 

(Xs are increasing.) 


The other two formats are wrap-around clockwise and wrap around 
counterclockwise, where the coordinates begin at the trailing 
edge, wrap around the leading edge, and end at the trailing 
edge* The clockwise case begins with the lower surface, and the 
counterclockwise case begins with the upper surface. The format 
shown below is essentially the same for both cases. NPTS is the 
total number of points in the dataset. 


FT! 


6 



PROFILE User's Guide 


LEGEND 

NPTS 

X Y 

X Y 


< CHARACTER* 80 > 

< 14 > 

< 2F10 . 0 > 

< 2F10. 0 > 

(May be X/C, Y/C; 

Xs are decreasing 
until the leading 
edge, then increasing) 


NOTES : 

* With PROFILE format the lower surface coordinates are 
optional, but a zero must be read for NL if no lower 
surface is included. 


* If both surfaces are present (PROFILE format), they must 
have the same leading edge point . The trailing edge points 
may differ. 

* The descriptive first line is termed a "legend" because it 
appears as a legend, not a title, on the airfoil plot that 
is obtainable with PROFILE. 


* There may be more than one set of coordinates in the same 
file (all formats). 

* Multi-element airfoils may be in any of the three formats, 
with a set of coordinates for each element. Unless each 
element is to be represented separately in the legend, one 
legend entry should be included with the first set of 
coordinates and a blank line should precede any following 
sets. 


* The 2F10.0 is not strictly necessary, as PROFILE uses 
list-directed reads rather than formatted reads. This is 
partly because the subroutine PRREAD is used for reading 
more than just airfoil coordinates. But 2F10.0 is a 
popular standard. 
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2.2 Control Inputs 


Control inputs common to all modes are represented by 
keywords and values in the input control file. PROFILE prompts 
for the name of this control file. A description follows, with 
the default value for each keyword in parentheses. Note that 
synonyms are permitted in a couple of cases. Note also that 
several keyword /value pairs may appear on one line, except for 
the multi-valued keywords LINE and NOFILE, which must not appear 
with other keywords. 


KEYWORD VALUES DESCRIPTION 

FORMAT (SAME) One of three formats for 

PROFILE or STANDARD output profile(s). May be in 
CLOCKWISE standard PROFILE format, 
COUNTERCLOCKWISE clockwise wraparound format 

or counterclockwise wrap- 
around. SAME means the same 
format as the input profile. 


LEADING_EDGE (BLUNT) 

SHARP 


LINE (SOLID) 

(DASH) 

DOT 

CHA INDASH 
CHAINDOT 
THICK 
SYMBOLS 


Type of leading edge. Para- 
metric splines are used for 
airfoils with blunt leading 
edges and ordinary natural 
splines are used for air- 
foils with sharp leading 
edges . 

Controls line types of 
curves on profile plots. One 
value may be included for 
each curve on the plot . The 
default is a solid line for 
the first curve, a dashed 
line for the second curve, 
then lines with symbols for 
any following curves. The 
first curve would typically 
represent the original pro- 
file, with the second curve 
representing the revised 
profile. Overriding the de- 
fault might be necessary 
when plotting multi-element 
airfoils, or when overlaying 
more than 2 airfoils. At 
most 20 curves are provided 
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for. Note: The keyword 

SYMBOLS refers to symbols 
without a line connecting 
them. 


MAXCURVATURE (+5.) Cutoff values for plotting 
MINCURVATURE (-5.) curvatures. Practice shows 

that +/-5. give useful plot 
scaling by ignoring the high 
curvature values near the 
leading edge. On the other 
hand, the user may well want 
to focus on the leading edge 
region . 


NOFILE 


PLOT 


THREED 


(NONE) 

DAT 

PLT 

TAB 

CRV 

YPP 

(BOTH) 

ORIGINAL 

REVISED 


(FALSE) or (NO) 
TRUE or YES 


Used to suppress any com- 
bination of the five output 
files generated by PROFILE. 
The values correspond to the 
file name extensions - see 
BASIC OUTPUT DESCRIPTION. 

Controls plotting of ori- 
ginal and revised profiles. 
The default is to plot both 
original and revised (if one 
exists) . 

For plotting of multiple 
stations from a 3D wing . 
The default is the 2D case, 
giving all airfoil curves on 
the same plot frame. When 
THREED is true, each airfoil 
found in the input dataset 
is plotted on a separate 
frame . 


XAXIS (6.4) 


XMIN 


(minima and 

XMAX 


maxima of 

YMIN 


the input 

YMAX 


coordinates) 


Length of x-axis in inches. 
The default is appropriate 
for an 8.5 x 11 page. 

Minimum abscissa of window 
Maximum abscissa of window 
Minimum ordinate of window 
Maximum ordinate of window 


< Any or all of the last four keywords may be used 
to plot a portion (window) of the profile plot. > 
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A sample input control file follows. Note that keywords 
vfliues mav be separated by blanks, commas, colons, equal 
Sgns Sr tlbs Y As LntfLed abovl keywords with bore than ore 
value should appear on a separate line. Keywords (and text type 
vSl"s) My be truncated to any set of unambiguous leading 

characters . 


LEADING_EDGE : SHARP FORMAT = STANDARD PLOT REVISED 

MINCURVATURE - -6. MAXCURV 6. _ 

MIS -16 3 XMIN - .2 XMAX: .8 YMIN -.03 YMAX ■ 

LINE = SOLID, SOLID, DOT, DOT 
NOFILE : DAT YPP 


THREED: NO 
.04 
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3.0 BASIC OUTPUT DESCRIPTION 


PROFILE produces up to five different output files. The 
NOFILE keyword in the input control file (previous section) can 
be used to suppress any unwanted files. All the output files 
are named with a common identifier chosen by the user, each with 
a different extension. The five standard extensions are listed 
below, using the default identifier PROFILE. 


PROFILE . DAT 
PROFILE . PLT 
PROFILE. TAB 
PROFILE. CRV 
PROFILE. YPP 


Revised airfoil data file (requested format). 
Plottable airfoil data file (QPLOT format). 
Tabulations and diagnostics, for printing. 
Curvature data file (QPLOT format). 

2nd derivatives file (PROFILE format). 


Each output file is described in more detail in the 
sections that follow. 


3.1 Revised Airfoil Coordinates 

The < identifier >. DAT file contains the airfoil coordinates 
that have been revised in some way. The file can be in any of 
the three formats described in the BASIC INPUT DESCRIPTION and 
can be used in subsequent runs of PROFILE. 


3.2 Plottable Airfoil Data 

The < identifier >. PLT file contains airfoil coordinates in 
the plottable form handled by QPLOT (Appendix B) . Normally, one 
plot frame is produced (but see THREED usage below). A plot of 
one single-element airfoil may contain up to two curves, one 
representing the original input coordinates, and another 
overlaid curve to represent the revised coordinates. The 
keyword PLOT in the input control file specifies which of the 
curves are plotted (i.e. , original, revised, or both). For 
modes 1 through 3, only the revised profile is plotted, since 
there should be no visual difference between the original and 
revised profiles in these cases. 

The plot may contain more than 2 curves (up to 20) in the 
case of a multi-element airfoil, or when several airfoils are to 
be superimposed for comparison. If multiple sections of a wing 
are "indicated via the THREED keyword, these are plotted 
separately on multiple frames. 
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A prompt is issued for a title which will be centered at 
the top of the plot. This title will also appear as the legend 
entry in the revised airfoil dataset (if one is requested), and 
is used to identify any revised data in the < identifier > .TAB and 
< identifier > .CRV files. An optional subtitle will also be 
centered under the title. 

A variety of line types may be specified in the input 
control file. (See BASIC INPUT DESCRIPTION for details on line 
types , particularly for the multi-element case.) A legend 
centered above the axes identifies the curve(s). 


3.3 Tabulations 

The < identifier > .TAB file contains tabulations of the 
profile coordinates as well as the values of 1st and 2nd 
derivatives and curvature at each point. When the profile has 
not been revised, tabulations for only the original profile are 
printed otherwise tabulations for both the original and revised 
profile are included. Details of the iterative computations are 
also printed here if REFINE or OPTIMIZE mode is being used. 
Other diagnostics written to the .TAB file include the number of 
profiles identified during the preliminary scan of the .DAT file 
(for overall data range purposes), and a record of the shape 
functions selected if the MODIFY option is used. 


3.4 Curvature Data 


The < identifier > .CRV file contains abscissas and curvature 
values in the form expected by QPLOT. Upper and lower surface 
curvatures will be plotted on separate frames, with up to three 
curves on a frame. Original curvatures are included for mode 0, 
original and revised curvatures for modes 1 through 5, or 
original, revised and target curvatures for mode 6. (For 

details on target ourvature distributions, see the OPTIMIZE mode 
in Chapter 4 . ) By default , the curvature values plotted will 
range between -5. and +5., but these limits can be redefined 
using the input keywords MINCURVATURE and MAXCURVATURE. (All 
values are written to the .CRV file; it is QPLOT 

suppresses those outside the chosen range during plotting. ) 
curvature data can also be edited and reused as input 
running PROFILE in the OPTIMIZE mode . 


that 

The 

when 
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3.5 Second Derivatives 

The < identifier > .YPP file is a table of 2nd derivatives at 
each point on the airfoil in standard PROFILE format. When the 
profile has not been revised (M0DE=0), derivatives for the 
original profile will be in the file, otherwise only derivatives 
for the revised profile are included. This avoids mistakes when 
the file is reused (probably after editing) as an input file for 
PROFILE'S REFINE option. In this mode, the Y" values are used 
to act as constraints on the curvature. See section 4.6 for 
more on the reuse of 2nd derivative values. 
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4.0 MODES OF OPERATION 


Although there are several inodes of operation available in 
PROFILE, all of them require the same basic inputs (i-©- * * n 
airfoil dataset and a control input file), and nearly all of the 
output options are available for each mode. The principa 
difference between modes is in the method(s) used to manipulate 
the airfoil geometry, which are detailed in this chapter. 


PROFILE can be run in only one mode at a time, selected 
from a menu. It may be necessary to run the program several 
times for one airfoil, depending on the options desired. These 
options are described in roughly the order of complexity. 


4.1 MODE=0 : Tabulate and/or Plot (No Modifications) 

This mode does nothing to modify the input coordinates 
(with the possible exception of writing the coordinates to disk 
in a format different from the input format). Select MODE-O 
simply to plot an airfoil and/or its curvature distributions, or 
to generate the 2nd derivatives needed for reuse (after editing) 
in a REFINE mode run. 


4.2 M0DE=1 : Rectify Leading Edge Definition 

This mode is used to "rectify" the leading edge definition 
in cases where the current leading edge is not the most forward 
point on the airfoil. For example, some NACA airfoils have 
negative abscissas on the upper surface and leading edge defined 
as (0 0 ) This causes the upper surface to be multi-valued, 
creating' problems with spline fitting and the addition of shape 
functions M0DE=1 finds the minimum abscissa on the airfoil and 
uses this as the leading edge point common to both surfaces, 
shifting the upper and lower surface points accordingly. The 
result is an airfoil with a "true" leading edge and 
monotonically increasing abscissas. 


4.3 MODE-2: Normalize/Denormalize 

This mode renormalizes airfoil coordinates. PROFILE 
determines the chord of the (first) input airfoil and prompts 
the user for whether to normalize the coordinates by this chord 
(resulting in an airfoil with a chord of 1 . ) or by some other 
value. If some other positive value is entered, a prompt will 
be issued for the coordinates of the leading edge point as well. 
If the chord value entered is negative, the (presumably 


! I 
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normalized) profile will be denormalized by the absolute value 
of that chord. 

Care must be taken when more than one airfoil is processed 
at a time, as with a 3D case or a multi-element airfoil. The 
same renormalization is applied to all airfoils in the input 
dataset . 


4.4 


M0DE=3: Redistribute Airfoil Coordinates 


This mode redistributes the airfoil abscissas in one of 
several ways, according to the option chosen from another menu. 
Either parametric or ordinary natural splines are used to derive 
the ordinates, depending on the value of keyword input 
LEADING_EDGE . The menu selections are shown below. 


MENU SELECTION 
-1 


0 , 1,2 


0 

1 

2 

3 


DESCRIPTION 

Here, the new abscissas are to be read from a 
file prepared by the user in standard PROFILE 
format. Ordinates may or may not be included, 
although they will be ignored. A line of text 
(even if blank) should be included at the 
beginning of the file in place of the legend 
entry. A prompt will be issued for the name 
of this file. 

These menu selections use the subroutine 
XGRID to redistribute abscissas within the 
range of the original coordinates. Prompts 
are issued for the desired number of ab- 
scissas to be redistributed on the upper and 
lower surfaces. Descriptions of each of these 
options follow. 

Uniform spacing. 

Sinusoidal bunching near the leading edge . 

Sinusoidal bunching near both the leading and 
trailing edges. 

This menu selection is reserved for some 
other method of redistribution in a user- 
supplied version of XGRID (in which case the 
program would require relinking). 


15 



PROFILE User ' s Guide 


4.5 MODE-4: Modify Airfoil Using Perturbing Functions 

This mode allows the user to perturb each of the airfoil 
surfaces by applying shape functions interactively, according to 
the menu below. 


SCALE 1 
TRAILING 2 
DROOP 3 
LEADING 4 
EXPONENTIAL 5 
SINE 6 
WAGNER 7 


Y = Y * PI 

Y - X ** PI 

Y = e ** ( -PI * X ) 

Y-(l. - X ) * * PI 

Y = X * * PI * ( 1. - X) * e ** C -P2 * X ) 

Y = SIN ** P2 ( PI * X ** (LOG( . 5 ) /LOG( PI ) ) ) 
Wagner function; PI - Order of requested term 


Where appropriate, prompts are issued for the parameters PI 
and P2 and for a multiplier of each function selected. Each 
surface is perturbed independently with as many shape functions 
as desired. There is an option to start over in case of error, 
before moving on to the next surface. The selected shape 
functions , parameters , and multipliers are echoed in the 
< identifier > .TAB file, which may be printed for future 
reference. See Appendix C for illustrations of the shape 
functions available. 
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4.6 MODE =5 : Refine Airfoil by Manipulating Y" 

4.6.1 Motivation and Analysis: 

The REFINE option that is now part of PROFILE arose from a 
requirement to adjust the thickness ratio of an airfoil that had 
been produced by the design-by-optimization program FL06QNM 14 
Simple scaling of all of the ordinates is known to cause loss of 
some of the properties that have been laboriously achieved 
through the optimization process. The curvature m the leading 
edge region, and perhaps the trailing edge region, plays a 
significant role in the flow characteristics, so a technique was 
sought which would retain the original curvature distribution as 
much as possible while producing a desired thickness. 

Curvature, K, is defined in terms of first and second 
derivatives, Y' and Y" , as follows: 

Y" . 


2 3/2 

( 1 + Y' ) 

The motivation for the method implemented as REFINE is that, if 
Y' is small, then curvature K is very like second 
Now Y" can be represented as a finite difference 
which is a linear combination of ordinates Y: 

2 ( DX Y - ( DX + DX )Y + DX Y ) 
i i-l i i-1 i i-1 i+1 

Y" = 

i 

DX DX ( DX + DX ) 
i i-1 i i-1 


where DX X - X 

i i+1 i 

Reading Eq. 4.6.2 the other way around suggests 
simultaneous equations in the UNKNOWN ordinates 
DESIRED values of the 2nd derivative (and hence, approximately, 
the curvature) on the right hand sides. Furthermore, airfoil 
thickness can be modified by some sort of scaling of ^ ® 
original ordinates. This suggests an additional set of 
equations. Thus an overdetermined system of linear equations 
for the ordinates can be set up, the optimal solution of which 
tends to satisfy the combined requirements involving thickness 
and curvature . 


expression 


(4.6.2) 


a set of 
Y(i) with 
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The basic idea of posing the problem as a linear least 
squares one proves to be more powerful than was initially 
envisioned, as will be explained below. But first, a look at 
the system that is set up and solved for N points defining one 
surface : 


Z(i) - scale(i) * Y(i) (i=2:N-l) (4.6.3) 


Z " ( i ) * Y" (i) (i=2:N-l) (4.6.4) 


where Z represents the revised form of ordinates Y; 

and scale is the (normally nonlinear) distribution of 

scale factors which yields the specified 
thickness ratio in an iteration to be de- 
scribed below. 


Note that end points i=l and i=N are omitted because the 
central difference formula does not apply there. This is not a 
problem, since the values for Z(l) and Z(N) must equal 
scaled )*Y(1) and scale(N) * Y(N) - not just approximately - in 
order that the two surfaces retain their common leading edge 
(and possibly trailing edge) point. 

The system shown as Eq. 4.6.5 is of order 2(N-2) by N-2, 
and is very sparse since the upper half is strictly diagonal and 
the lower half is tridiagonal. The two terms "B" in the right 
hand side represent the desired values of Y" for i=2 and i=N-I , 
adjusted in the usual way by the appropriate multiples of the 
known values for Y at i=l and i=N. SY represents scaled) *Y(i) . 
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i=2 


i=N-2 

i=2 


i=N-2 


> 

’ X 


■ Z ■ 


SY ‘ 


X 


z 


SY 


X 


z 


SY 


X 


z 


SY 

> 

X 


. z . 

« 

SY 

> 

X X 


B 


XXX 


Y" 


XXX 


Y" 


XXX 


Y" 

— > 

. X X - 


. B . 


(4.6.5) 


4.6.2 Weighting the Y" Equations: 

Since values of Y" tend to be considerably larger than 
values of Y (particularly near the leading and trailing edges), 
the two halves of Eq. 4.6.5 require equilibrating in some sense. 
This is achieved by applying multiplicative weights less than 
unity to the lower half, so that Eq. 4.6.4 becomes 


weight(i)*Z" (i) = weight (i) *Y" (i) (i=2:N-l) (4.6.6) 


Furthermore, nonuniform weighting of the Y" values is advisable, 
just as nonuniform scaling of the Y values is necessary, but for 
different reasons. In Eq. 4.6.3, scale(i) must vary from a 
value of 1.0 at leading and trailing edges to that factor at the 
point of maximum thickness which gives the desired thickness. 
In the case of the second derivative weights, smooth variation 
over an order of magnitude has been found to give best results • 
typically 0.004 at leading and trailing edges up to 0.04 at some 
interior point (yet to be specified) . (This sort of 
nonlinearity reflects to some extent the fact that Y" is much 
smaller in magnitude in the mid-section than it is fore and aft 
more than an order of magnitude, typically. The goal is to 
retain the curvature as much as possible towards the extremes, 
so the -outer equations of Eq. 4.6.6 still tend to dominate in 
terms of residual over the inner equations. ) The weighting 
prevents the residuals of Eq. 4.6.3 from being swamped by those 
of Eq. 4.6.4. 
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4.6.3 Smoothing the Y" Right -Hand- Side Values: 

Selection of the point with maximum weight cannot be 
described until an unanticipated benefit of this computation is 
explained. The reader may be surprised if he plots a few 
tvpical curvature distributions: they can be remarkably spiky 
and far from smooth. It makes little sense to strive to retain 
these spiky distributions. Smoothing out the Y" values used on 
the right-hand-side of Eq. 4.6.6 is highly desirable. And here 
is where the aptness of the term REFINE may be appreciated: we 
have a means of tailoring the Y n distribution (and hence the 
curvature) to any reasonable, smooth, shape we choose. Since 
smoothness and general shape are more fundamental than absolute 
values here, it does not matter that all equations in the 
overdetermined system cannot be solved exactly. Moreover ( and 
this answers any objection that Y' is NOT negligible in the fore 
and aft regions of most interest), it does not matter that we 
are working with Y” rather than curvature K. For if the 
right-hand-side values used in Eq. 4.6.6 are related to, or 
varied systematically from, the original Y values, and t e 
revised values Z for ordinates Y bear some relationship to the 
original Y values (as they do by construction), then the first 
derivatives of the Z values, and hence the revised curvature 
values, cannot help but bear some relationship to the original 
values - smoothed out in the qualitative way that is now seen to 
be desirable . 


4.6.4 Nonlinear Weights and Scales: 

Returning to the question of where to choose the peak 
weight for the Y" values, the answer is that point where some 
spike or unwanted bump in the original curvature distribution 
showed up. We have a means of zeroing in on the offending 
features of the distribution (not necessarily all in the same 
run, though!). 

Thus the kind of nonlinearity sought (for both the weights 
and the scales) involves known values at the end points and a 
known value at a known interior point - the center of some 
glitch in a region of the curvature in the case of the weights, 
or the abscissa where the airfoil thickness is greatest in the 
case of the scale factors. A function with the necessary 

properties is given by 


P3 * ( SIN ( PI * X ** PI ) ) 


* * 


P2 


(4.6.7) 


where PI - LOG(0.5) / LOG (X-at -peak ) , which serves to locate 
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the peak of the sine in the desired place; 

P2 - a "width" parameter, where 2.0, 3.0, or 3.5, say. 

give more and more "peaky" distributions, while 
1.0, 0.5, ... give broader distributions; and 

P3 « a multiplier to alter the sine's peak value, which 
is 1.0 for any PI and P2. 


The scale factors of Eq. 4.6.3 require further explanation. 
If the equation for the point where maximum thickness occurs 
were solved exactly, it can be shown that the desired thickness 
would be obtained by using 

P3 = i.o - <desired thickness > / ‘original thickness > 

and also using 1.0 - P3 * sin . . . rather than P3 * sin ... 
Note that this scaling is applied to either surface or to both, 
since S*( YU - YL ) = S*YU - S*YL in the case of both surfaces. 

However, none of the equations of the overdetermined system 
is solved exactly. This means an iteration is necessary, 
whereby the < desired thickness > above is treated as a target 
which is adjusted away from the true desired thickness until a 
match is obtained to within some suitably small tolerance. The 
algorithm for adjusting the "target" involves two "calibration" 
computations, after which linear interpolation between the two 
most recent estimates suffices till convergence. An oscillation 
may be observed between two abscissas showing maximum thickness 
if the thickness calculation is confined to the original, 
discrete, abscissas. This tends to slow the convergence, but 
should not be a problem. (More continuous computation of 
thickness using a spline on one surface must be used if the two 
surfaces are not defined at the same abscissas. However, the 
present implementation switches to the more correct calculation 
only if it detects different sets of abscissas.) 


4.6.5 Setting Up the Target Y" Values: 

Preparation of the target Y" values also requires 
explanation. This has been made as convenient as possible. If 
the curvature distributions are already smooth and about as good 
as the user knows how to ask for, then REFINE can compute and 
use these with no other inputs. 

Alternatively, some single region at a time may be targeted 
for smoothing simply by responding to prompts for a constant 
value of Y" over some Z-range. This works well for spikes, or 
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broader- bumps, in the mid-section curvature, and can also be 
handy for constraining curvature away from zero, where the 
airfoil may be excessively flat or reflexed. It can also work 
in noisy regions near the leading or trailing edges, since 
constant Y" (which tends to lead to smooth Y and hence Y') can 
still lead to smoothly varying curvature as required. To repeat 
an earlier statement : it does not fiiatter that we are 
constraining Y" rather than curvature. 

More generally, REFINE has the option of reading a table of 
Y" values so that more than one unsatisfactory region can be 
refined in a single run. This table is prepared in exactly the 
same format as for the airfoil itself, and may start out as the 
table saved by an earlier run of PROFILE. The table does not 
have to be complete: linear interpolation is done for abscissas 
in the table range (meaning chunks can be edited out where 
linear variation of Y" is considered acceptable), while the 
original Y" values are used for abscissas outside the table 
range. This kind of tailoring of the second derivative targets 
may take practice to get a feel for what is achievable and what 
is not. It is certainly flexible, but the user should keep in 
mind the simpler mechanism for "zapping" isolated spikes or 
bumps in the curvature described previously. 


4.6.6 Solution of the Overdetermined System: 

While the system set up and solved can be thought of as two 
groups of constraints - one on ordinates, one on second 
derivatives, it is treated here as an ordinary overdetermined 
system of equations with no constraints in the linear 
programming sense. In the original implementation, no attempt 
was made to take advantage of the sparsity of Eq. 4.6.5. The 
simplicity of a good direct method for dense systems, using the 
QR factorization [6] as in subroutine HDESOL, was preferred to 
existing alternatives for sparse systems that were likely to be 
iterative . 

However, this meant a lot of storage for the 2M * M matrix 
(M=N-2), and a lot of computation for typical values of M - a 
computation that is usually iterated 3—5 times till the desired 
thickness is achieved. (The storage requirement was about 300 * 
150 = 45000 words for up to 150 points per surface to be 
handled while the CPU time per iteration of the refinement 
calculation (both surfaces refined) was measured at about 16 
seconds on a VAX-1 1/780 for the typical value of 72 points per 
surface . ) 
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A specialized scheme has therefore been implemented for 
this diagonal -plus-tridiagonal structure. Subroutine DTDLSQ 
still forms the QR factorization, but much more efficiently than 
the general-purpose HDESOL . (The coding for use of HDESOL has 
been left in the source for REFINE, commented out, as an aid to 
understanding the set-up of the overdetermined system.) The 
above storage requirement has been reduced to 4 * 150 = 600 
words, while the CPU time per refinement iteration has become 
about 0.8 seconds on the VAX (both surfaces, 72 points per 
surface) . 

Note that working with the QR factorization [6] of this 
matrix A is certainly preferable to setting up and solving the 
too-often-seen Normal Equations (involving A-transpose * A). 
Not only would the latter method require double the precision 
(word-length) to give an accuracy comparable to that of the 
single-precision computation used here. 


4.6.7 Some Further Notes on Usage: 


The available bef ore-and-af ter plots of curvature, plus the 
tabulations, should help the user to decide how to vary some of 
the interactively-entered inputs. It may take several tries to 
achieve acceptable results. Sometimes one should use REFINE on 
the results of the last REFINE application; in other cases, it 
may be better to throw away a result and back up for a retry. 


Another point to remember is that the default range for the 
curvature plots is intended for a good presentation of the mid- 
section by suppressing the handful of large values that appear 
near the leading edge and possibly trailing edge. The user may 
not want to ignore these larger values - plotting them by 
expanding the range defined by MINCURVATURE and MAXCURVATURE is 
advisable, along with observing them in the tabulation. 


Note that finite differencing is used for estimating Y' and 
Y" , and hence curvature, rather than the analytic derivatives 
that splines can provide. The experience of the authors is that 
interpolating splines tend to produce surprisingly wavy 

distributions for the higher derivatives, even for "good" 
airfoils. The sign of Y" can be. very sensitive to small 
changes. Finite differences, on the other hand, exhibit 
steadier behavior. More importantly, since the essence of the 
method is to use a central difference expression for Y" to 
define the lower half of the left-hand-side matrix of Eq. 4.6.5, 
consistency demands that the right-hand-side terms also be 

derived using central differences. Thus PROFILE eschews splines 

throughout for tabulation purposes, and for computations, though 
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it does use them for redistributing data points. 

If the airfoil is defined on too fine a mesh, using the 
traditiona? 2F10.6 format, zeros can easily show up for the Y" 
estimates near the leading edge where they are clearly no 
reasonable. This is due to catastrophic cancellation in tie 
presence of too few significant digits. Going to, say, E15.6, 
helps superficially, but makes little sense because airfoils 
cannot be constructed to one part in a million. Tll f moral. 
Avoid excessively fine airfoil meshes: 192 points in the range 

[0. , 1 . ] is too many. 


4.6.8 Summary of Interactive Inputs Used in REFINE Mode: 

Short descriptions of the quantities for which REFINE 
prompts the user are given below. All have defaults where 
possible, as indicated by <CR> (Carriage Return). The user is 
urged to experiment with various combinations till some sort of 
feel for cause and effect is attained. Only then can the power 
of REFINE be exploited effectively. 


* Desired % thickness: <CR> retains present thickness. 

* Width param. for Y: Affects the nonuniform scaling 

applied to the ordinates (both 
surfaces . The default is 2.0. Larger (3.0 — 4.0) 
tends to retain leading/trailing edge shape more 
while 1.0 would constrain fore and aft less. 

* Input Y" table: <CR> means there is none, else 

the file name is entered. This 
file should be in the same "PROFILE" format as 
the airfoil coordinates. It can cover any range 
of abscissas. (Linear 'interpolation is used.) It 
may be an edited version of the file from a pre- 
vious run of PROFILE, or it may be much cruder. 

WARNING: If you delete points as you smooth out 
this distribution, it is easy to forget to 
update the integer number of points at the 
beginning of the data for each surface. Incor- 
rect counts here can lead to obscure errors . 

* Constant Y" value: <CR> means no such constraint; 

retain existing curvature val- 
ues as much as possible. Otherwise, a value of 
Y" entered will be sought in the abscissa range 


! Ii 
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that is prompted for next. 

* Corresp. X range: Enter low and high X values on 

the same line. Allow for the 
fact that strict inequalities are used when the 
program tests for being within this range. E.g. « 
Enter .39 .61 or .39,. 61 if you intend the 
constraint to apply over [0.4, 0.6]. 

* Width param. for Y" : Default is 3.0. Affects non- 

uniform weighting of the equa- 
tions representing 2nd derivative constraints in 
the overdetermined system being solved. Since 
the actual values of the 2nd derivatives being 
sought also act in a weighting sense, effects of 
this variable are not easy to predict. Values of 
2.0 or 1.0 should tend to let Y" change more. 

* X for peak Y" weight: Default is 0.5. The absolute 

* values of Y" are so much big- 

ger than those of Y that they all need to be 
scaled down in the system being solved. If you 

are trying to flatten the curvature plot in some 
region, then pick the center of the region for 
this input. Otherwise, use the mid-chord value. 

* Y" weights, X/OO.l: Default is 0.004. See next. 

* peak Y" weight: Default is 0.04. These provide 

for the fact that the absolute 
values of Y" are typically smaller in the mid- 
section than near the leading /trailing edges, so 
they should be weighted more, especially in view 
of the fact that any Y" constraints applied are 
typically in the mid-section. See above. 
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4.7 M0DE=6: Optimize Airfoil by Manipulating Curvature 


4.7.1 Motivation and Analysis: 


The OPTIMIZE option incorporated in PROFILE is a natural 
generalization of the ideas embodied in the REFINE option. This 
mode of perturbing an airfoil works directly with the curvature 
distribution rather than the second derivatives. The price 
paid is that linear techniques can no longer be employed. 
Nonlinear optimization is typically much more expensive in terms 
of CPU time. On the other hand, OPTIMIZE mode offers a design 
method that is much cheaper than working with pressure distribu- 
tions which involve at least one flow field solution for each 
evaluation of the objective function being minimized. The hope 
is that at least some of the work it takes to optimize airfoils 
can be done on a modest scale by tailoring the curvature distri- 
bution first. 


The implementation is analogous to that of optimizing flow 
codes such as FL06QNM * : set up some target distribution and 

seek to match it as well as possible in the least squares sense 
bv applying appropriate shape functions or "bumps". As is true 
with flow optimization, constraints of some kind are almost es 
sential. Here, it is most likely to be the maximum thickness 
that has to be constrained somehow. Given that the minimiz- 
ation algorithm used, QNMDIF* , is an unconstrained one, we have 
two choices : use a penalty function to keep thickness within 

reasonable limits (and then perhaps use REFINE mode on the 
result to give the preferred thickness precisely): or hope for 

the best by solving without constraints and, again, relying on 
REFINE to adjust the result. If the target curvatures are at 
all reasonable, the second approach may well be adequate. 


Note that there is no "multiple design point" complication 
(referring to the desirability of optimizing airflow character 
istics at more than one flight condition at the same time). In 
fact no attempt is made to manipulate more than one surface a 
a time - the upper and lower surfaces can be optimized mdepend 
ently in separate runs of PROFILE, with little loss of general- 
ity and significant simplification. 


Prospective users of the OPTIMIZE option should be aware 
that it will probably take a lot of tries with different shape 


* FL06QNM [4] links a 2D potential flow code, FL06, with a 
reliable optimization package known as QNMDIF [1] plus a 
scheme for estimation of optimal finite differencing 
intervals. FL06QNM was developed under the same contract as 
PROFILE . 


I 
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functions, initial airfoils, target curvature distributions, and 
scale factors, before an acceptable result is achieved from a 
sufficiently well-posed, well-scaled optimization problem, . The 
initial airfoil is included here as a variable, because if it is 
not smooth to start with, perturbing it with smooth shape 
functions cannot be expected to produce smooth curvature distri- 
butions close to the (presumably smooth) target distribution. 
Use of the REFINE option first may be desirable, which raises 
the question of staying with the reliability and greater simpli- 
city of REFINE in favor of OPTIMIZE' s uncertain behavior. (The 
optimizing package may not converge to a minimum of the object- 
ive function, though it will almost certainly reduce its value 
below the initial value. But a lowered objective function does 
not necessarily represent an acceptably smooth solution, as in 
the case of oscillations observed occasionally when just Wagner 
functions are applied.) 

To this, the response must be that one does not know until 
one tries both possibilities. The OPTIMIZE option is presented 
as a potentially powerful tool, exercised by the authors on only 
a limited number of applications at the time of writing. Feed- 
back from adventurous users is welcomed. 

An outline of the OPTIMIZE method used in PROFILE follows: 

* Select upper or lower surface, and save the original 
coordinates for repeated application of varying shape 
functions . 

* Display present thickness ratio; prompt for desired one. 
Defaulting means either to retain the same thickness, or 
that the user doesn't care; in the latter case, the 
prompt for the penalty parameter to be used to impose a 
thickness constraint should be defaulted (to zero, mean- 
ing no such constraint). 

* Set up the user-defined "bumps" , either by parsing a 
file of keywords and values prepared ahead of time (more 
on this below), or by choosing the first N Wagner func- 
tions, which are easily definable without such a file. 

* Read the target curvature distribution, which is most 
likely an edited form of an output file from an earlier 
run of PROFILE. 

* Extract the active variables from the given set - not 
all of the parameters involved will be free to change 
some may be specified as inactive or fixed. 
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* Compute the initial value of the objective function, 
and estimate optimal finite differencing intervals. 

* Minimize the objective function (quasi-Newton method, 
without analytic derivatives). 

* Repeat the best function evaluation, and update the 
airfoil surface permanently. 

* Display the thickness actually achieved; perform the 
rest of PROFILE'S plotting and tabulations. 


4.7.2 Preparing the Target Curvature Distribution- 

One of PROFILE'S options is to generate a plottable file of 
curvature data. Typically, the plot for one of the surfaces, or 
both will exhibit imperfections that the designer would like to 
eliminate. PROFILE offers two choices: the REFINE option, which 
which uses second derivative distributions rather than curvature 
and is described elsewhere; and the OPTIMIZE option which works 
directly with curvature. The designer edits the imperfect values 
to smooth them out in some way, assisted by the original plot, 
then reruns PROFILE in OPTIMIZE mode, hoping that the new plots 
of curvature (showing revised and original distributions super 
imposed) will show a good match. 

Remember that OPTIMIZE mode can work with only one surface 
at a time. The input file of target curvatures is assumed to 
contain only the relevant surface. It is also assumed to retain 
the title and two axis labels that the original QPLOTable file 
contained. The title is transmitted to the output tabulation, 
along with the target curvature values; the labels are ignored. 
(More precisely: the target file may contain both surfaces, but 
only the first surface is read, so it must be the appropriate 

one . ) 


There is a question as to how many of the chordwise values 
should be retained in the target curvature distribution. There 
are considerations such as the fact that the values n ® ar t ^ e 
leading edge, being larger in magnitude, may tend to dominate 
the sum-of-squares objective function. (Extreme leading an 
trailing edge values, however, are suppressed by PROFILE in the 
plottable files and tabulations.) In the absence of a s °keme 
for weighting the elements of the sum of squares (which could be 
implemented in much the same way has it has been in the REFINE 
option) it is probably best to retain target values for all 
interior abscissas, if only to avoid surprises in any region 
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omitted . 


4.7.3 Selecting the Shape Functions: 

Appropriate selection of the shape functions is probably 
the most difficult aspect of using PROFILE'S OPTIMIZE option. 
There is no single choice for any given situation. Experience 
is desirable to help pick a subset that should do the job. In 
general, all of the bumps selected should be distinctly differ- 
ent in some sense - it is unwise to include two functions that 
have similar effects, as this adds unnecessary variables and 
blurs the minimum from the optimizing algorithm's point of view. 
If the changes being sought in the curvature are localized in 
some way, then some of the functions can be localized or center- 
ed similarly. 

PROFILE offers two basic sets of perturbing functions: six 
"bump" functions [3], and a family known as the Wagner functions 
[5]. All of these are illustrated in Appendix C. The 
perturbing functions to be used may be specified in two ways: 
in a previously-prepared file containing keywords and values; 
and a simpler way that is convenient if only the Wagner func- 
tions are to be used. (See below for more on the use of the 
Wagner functions.) 

Specification of the bump functions to be used has been 
made as mnemonic and straightforward as possible. There is no 
need for hard-coding the set of bumps - a keyword-driven input 
scheme is used, involving preparation of a small file that will 
look something like the following sample: 


BUMP: SINE 

CENTER: 0.5 STATUS: INACTIVE 

WIDTH: 3.0 STATUS: INACTIVE 

MULTIPLIER: 0. STATUS: ACTIVE SCALE: 100. 


BUMP - exp 
POWER = 15. 
WIDTH - 10. 
MULT = .001 


STATUS - fixed 

STATUS - fixed SCALE = 1. 

STATUS « variable SCALE: 10. 


29 



PROFILE User ' s Guide 


Points to note: 


* Format is free, with several possible delimiters; blank 
lines are optional, but good for separating bumps. 

* Keywords need only be long enough to be unambiguous. 

* Either upper or lower case is permitted. 

* Keywords for a given bump must begin with the BUMP key- 
word, which must be the only keyword on the line. 

* The ordering of subsequent lines describing that bump ' s 
variables or parameters is unimportant , but if some are 
omitted, they will be detected as undefined, and execu- 
tion will halt. There is NO ATTEMPT TO DEFAULT values. 

* Variable names for a given bump (e.g. WIDTH) must be 
the FIRST keyword on a line, one per line. 

* STATUS and SCALE may be in either order within a line . 

* SCALE is optional if STATUS is FIXED/ INACTIVE/ CONSTANT . 

* SCALE is defaulted if STATUS is ACTIVE/FREE/VARIABLE 
and no entry is given. (1.0 is the default.) 

* End-of-file signals end of data. Some special keyword 
would probably be necessary on a mainframe. 


The terminology requires further explanation. First, all 
of the perturbing functions (including the Wagner functions) 
have MULTIPLIERS which are almost certain to be included among 
the optimizing variables. Second, all of the functions are 
additive - multiples of the basic shapes are added to the 
original airfoil. Third, each bump type has a descriptive name, 
and is a function of a normalized abscissa x and of parameters 
with descriptive names, as shown below. (The SCALE function is 
an exception to all of these statements - it simply scales the 
ordinates y directly.) 
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Bump Name Parametric Expression 

DROOP exp ( -WIDTH * x ) 

EXPONENTIAL x**POWER * (1-x) * exp( -WIDTH*x) 

LEADING (1-x)** POWER 

TRAILING x * * POWER 

SINE ( sin ( pi * x**( log( 0 . 5) /lOg(CENTER) ) ) )**WIDTH 

SCALE Y * FACTOR 

WAGNER ( t + sin (t) )/pi - sin (0.5*t)**2 (N-l) 

( sin ( N* t )/N + sin ( (N-l)*t) )/pi (N > 1) 

where t = 2 * asin ( sqrt(x) ) 


4.7.4 Simplified Specification of Wagner Functions: 

If the only shape functions that the user wishes to apply 
for a particular run are the first N Wagner functions [5], then 
the keyword-driven file input description can be dispensed with 
since the set-up for these functions is easily generated. The 
only likely reason that the file input might still be preferred 
is that the given scaling of the variables may prove less than 
satisfactory. The present automated scheme generates all of the 
scale factors from a "seed" entered by the user. The relative 
scaling for N = 1,2, 3, 4, 5,... at the time of writing is 

1,1,1,1,1,3,3,3,3,3,10,10,10,10,10,30,30,30,30,30,100,100, 

based on very little empirical evidence. Such values should be 
established more systematically than experience has permitted so 
far, though optimal scaling in general is an open question. See 
"Choice of Scale Factors" below. 


4.7.5 Quasi-Newton Method For Nonlinear Optimization: 

The nonlinear optimization algorithm adopted here is that 
developed by Kennelly for program FL06QNM [4] , itself a refine- 
ment of the QNMDIF implementation of Gill, Murray, et al . [1]. 

Only a brief outline of this revised quasi-Newton method can be 
provided here. QNMDIF, in combination with the algorithm of 
subroutine FDSTEP (below), is a modern, robust, and efficient 
package for the unconstrained minimization of functions (not 
necessarily sums of squares) of N variables. It uses finite _ 
difference approximations for first derivatives, which in turn 
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are used to build, through successive updates, a matrix having 
Hessian-like properties. (The Hessian for afunotron of N 
variables is the NxN matrix of second derivatives . ) The 
approximate gradient and Hessian are used to compute a Nevton- 
S™ search direction at each step, with descent ensured by 
controlling the positive-definiteness of the matrix via its 
Cholesky (LDL') factors. 

A safeguarded line search routine chooses adequate step 
lengths efficiently, with account taken of functions that are 
calculable to less than full machine precision, and automatic 
switching between forward and central differencing in response 
to function behavior. 

This package is appropriate for the present application 
since analytic derivatives of the objective function are almost 
as out of the question as they are in the case of target press 
ure distributions. Furthermore, a good general-purpose package 
tends to behave as well as, or better than, typical packages 
specially adapted to sum-of-squares-type functions. Also, the 
relative simplicity of an unconstrained algorithm is preferred 
£o a general constrained algorithm, since the only likely con- 
straint is one on thickness, which can be handled with a penalty 
function (below). 

Experience so far suggests that single precision C32-bit) 
arithmetic is adequate for QNMDIF to solve reasonable OPTIMIZE 
problems. The significant QNMDIF parameter ' £oi 

no provision for interactive ohanges include ETA-0.1, TOL- • 001 , 
EPSMCH=5 . E-8 , and EPSOBJ=EPSMCH* 100 . The initial Hessian-like 
approximation is a diagonal matrix containing estimates for the 
subset of second derivative elements provide by the preliminary 
CENDIF/FDSTEP calculations (based on algorithm FD of [2]), 
are not optional either, as they play a large part ** ?™ DIF S 
success by estimating optimal finite differencing intervals. 

An arbitrary upper limit of 100 optimization iterations is 
imposed. This can lead to many hundreds of function evaluations 
if the problem is poorly posed, so the airfoil designer is urge 
to be careful in setting up the OPTIMIZE option. 

Restarting an OPTIMIZE solution is always a possibility. 
The improved finite difference intervals may make a difference 
on difficult problems. However, no provision has been made for 
saving the accumulated curvature information represented by the 
current Cholesky factors L and D, as this application of QN D 
(to a relatively well-behaved function) is most likely to find 
solutions in a single run in a modest amount of CPU time. 
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4.7.6 Choice of Scale Factors: 

The scale factors applied to active variables are intended 
to assist the optimizing algorithm, which behaves best if the 
solution lies in the unit sphere and the rows of the Hessian 
matrix of second derivatives at the solution are approximately 
equilibrated. (Ideally, the function should not exceed unity in 
the unit sphere, but it may not be possible to meet any of these 
criteria.) Many aspects of the scaling problem are discussed in 
[1] and [2]. To some extent, less-than-ideal scaling is handled 
satisfactorily by the dynamic selection of intervals for finite 
differencing by routines CENDIF and FDSTEP. 


4.7.7 Choice of Target Thickness And Penalty Parameter: 

OPTIMIZE prompts the user for a target maximum thickness 
ratio (as a percentage of chord). Unlike the REFINE option, it 
cannot guarantee this target thickness precisely. The target is 
used in a penalty function of the form 

< penalty parameter > * ( current % T/C - target % T/C )**2 

which is added to the rest of the sum-of-squares-type objective 
function. The user can seek the same thickness as the original 
airfoil by entering just a carriage-return (but it is unlikely 
to be retained precisely) . The user can also choose to suppress 
the constraint by entering either zero or a carriage-return in 
response to the prompt for the penalty parameter. This can save 
CPU time, because the thickness calculation is then avoided. In 
practice, however, the presence of a constraint of this nature 
can actually speed the convergence of the optimizing algorithm, 
which isn't to say that the target curvatures will be matched 
any better. QNMDIF should find a minimum, but there is nothing 
to stop the user from reformulating the problem if he doesn't 
like the best fit that was found. 

In principle, successive runs of OPTIMIZE with increasingly 
larger values for the penalty parameter can produce the desired 
thickness, but the availability of the much more precise, and 
cheap, REFINE option argues for an OPTIMIZE run with, say, 
penalty parameter = 100. , followed by a refinement run (with 
numerous options open) if the thickness is required to be some 
precise value. 
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4.7.8 Comparison of Curvature With The Target Distribution: 

The same table-look-up with linear interpolation that was 
implemented for the REFINE option's target y" requirements is 
used in routine CFDISTRIB . This means that, in principle the 
target and current curvature distributions do not have to be de 
fined at the same abscissas, though in practice it is most like 
lv that the two sets of abscissas will be identical. No attempt 
has been made to weight the elements of the sum of squares, in 
the interest of simplicity. 


4.7.9 Displaying Results: 


The user is prompted for minimal or comprehensive output of 
©NMDIF's progress. Normally, brief output should be all that is 
needed. Upon termination of PROFILE, the recommended sequence 
is as follows (where < ident > refers to the identifier prompted 
for at the beginning of the PROFILE run) : 


QPLOT < ident > . PLT 


revised-versus-original airfoil plot 


QPLOT < ident > . CRV ... upper/lower, old/new curvature plots 

PRINT < ident >. TAB . . . tabulations of the original airfoil, 

the target information, bump function 
description, CENDIF results, QNMDIF 
results, and the revised airfoil. 
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Fancy QPLOT Example 

Look, look! A rather long-ish subtitle is now optional. 

□ Observed (raw data) 

Macroeconomic Theory 
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All the somple results shown here contain tabulated output 
( i . e . <i dent i f i er>. TAB) files. Since these results include 
original and/or revised airfoil coordinates, derivatives and 
curvatures, the other output files are not shown in most cases. 
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Tabulated results: NACAMODE0 . TAB 
Curvature QPLOT file: NACAMOOE0 . CRV 
2nd derivatives file: NACAMODE0 . YPP 
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A normalized airfoil geometry input data set is assumed 
when running with MODE—3 or greater. The NACA 6-series airfoil 
of the previous example is normalized here. 
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Program PROFILE: REFINE mode used on Ames-01 Rotor Section 


ORIGINAL PAGE SS 
OF POOR QUALITY 
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PROFILE User’* Guide (Appendix D. Sample 5) 


Program PROFILE; REFINE mode used on Ames 01 Rotor Section 

Curvature smoothed via target y" distribution; new T/C=10% 

ORIGINAL AMES 01 
• is e d jpro fil e 



x/c 




upper surface curvature 


PROFILE User'* Guide (Appendix D» Sample 5) 


Program PROFILE: REFINE mode used on Ames-01 Rotor Section 

Curvature smoothed via target y" distribution; new T/C=10% 
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Program PROFILE: REFINE mode used on Ames-01 Rotor Section 

Curvature smoothed via target y" distribution; new T/C-10% 


□ = Original 
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The example of the OPTIMIZE option starts with the result 
from the previous exomple - a REFINEd Ames-01 rotor section. 
This may not be very representative in that the result achieved 
here, on something very like it, could have been achieved more 
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PROFILE User’s Guide (Appendix 0, Sompie 6) 


Program PROFILE: OPTIMIZE mode on REFINEd Ames-01 Section 

Target curvature distribution; no constraint on thickness 

Program PROFILE: REFINE mode used on Amea-Ol Rotor Section 
Revised 
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PROFILE User's Guide (Appendix D. Sample 6) 




Program PROFILE: OPTIMIZE mode on REFINEd Ames-01 Section 

Target curvature distribution; no constraint on thickness 

□ =Original 
G= Target 
A = Revised 
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PROFILE liter* s Guide (Appendix D, Sompie 6 ) 


Program PROFILE: OPTIMIZE mode on REFINEd Ames-01 Section 

Target curvature distribution; no constraint on thickness 
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